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Summary I propose a path integral description of the Su-Schrieffer-Heeger 
Hamiltonian, both in one and two dimensions, after mapping the real space 
model onto the time scale. While the lattice degrees of freedom are classical 
functions of time and are integrated out exactly, the electron particle paths are 
treated quantum mechanically. The method accounts for the variable range of 
the electronic hopping processes. The free energy of the system and its tem- 
perature derivatives are computed by summing at any T over the ensemble of 
relevant particle paths which mainly contribute to the total partition function. 
In the low T regime, the heat capacity over T ratio shows un upturn peculiar 
to a glass-like behavior. This feature is more sizeable in the square lattice than 
in the linear chain as the overall hopping potential contribution to the total 
action is larger in higher dimensionality The effects of the electron-phonon 
anharmonic interactions on the phonon subsystem are studied by the path 
integral cumulant expansion method. 

1 Introduction 

There has been a growing interest towards polarons over recent years also in 
view of the technological potential of polymers, organic molecules [IJ [H |3l 
[H O [7], carbon nanotubes [H] and high T c superconductors [3] in which 
polaronic properties, essentially dependent on the type and strength of the 
electron-phonon coupling jTUl [HI Q~2] , have been widely discussed. 

One dimensional (ID) systems with a half filled band undergo a structural 
distortion [13] which increases the elastic energy and opens a gap at the Fermi 
surface thus lowering the electronic energy. The competition between lattice 
and electronic subsystems stabilizes the ID structure which accordingly ac- 
quires semiconducting properties whereas the behavior of the 3D system would 
be metallic-like. 

Conjugated polymers, for example polyacetylene, show anisotropic electri- 
cal and optical properties [T3] due to intrinsic derealization of it electrons 
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along the chain of CH units. As the intrachain bonding between adjacent CH 
monomers is much stronger than the interchain coupling the lattice is quasi- 
1D. Hence, as a result of the Peierls instability, irans-polyacetylene shows an 
alternation of short and long neighboring carbon bonds, a dimerization, ac- 
companied by a two fold degenerate ground state energy. Charge injection in 
polymers induces a lattice distortion with the associated formation of local- 
ized excitations, polarons and/or charged solitons. While the latter can exist 
only in irans-polyacetylene, the former solutions are more general since they 
do not require such degeneracy. 

The Su-Schricffcr-Heeger (SSH) model Hamiltonian [15] has become a suc- 
cessful tool in polymer physics as it hosts the peculiar ground state excitations 
of the ID conjugated structure and it accounts for a broad range of polymer 
properties [16l \17\ [18] . In the weak coupling regime a continuum version [19] 
of the SSH model has been developed and polaronic solutions have been ob- 
tained analytically [5D] . Although the ID properties of the SSH model have 
been mainly investigated so far [2T] [22 [23] [IJ [13 [26] [27] [28] , extensions to two 
dimensions were considered in the late eighties in connection with the Fermi 
surface nesting effect on quasi 2D high T c superconductors [29]. Recent nu- 
merical analysis [30] have revealed the rich physics of 2D-SSH polarons whose 
mass seems to be larger than in the ID case, at least in the intermediate 
regime of the adiabatic parameter [31] . 

As a fundamental feature of the SSH Hamiltonian the electronic hopping 
integral linearly depends on the relative displacement between adjacent atomic 
sites thus leading to a nonlocal e-ph coupling (32] [33] with vertex function de- 
pending both on the electronic and the phononic wave vector. The latter 
property induces, in the Matsubara formalism [34], an electron hopping as- 
sociated with a time dependent lattice displacement. As a consequence time 
retarded electron-phonon interactions arise in the system yielding a source 
current which depends both on time and on the electron path coordinates. 
This causes large e-ph anharmonicities in the equilibrium thermodynamics of 
the SSH model [35]. Hopping of electrons from site to site accompanied by a 
coupling to the lattice vibration modes is a fundamental process determining 
the transport [36] and equilibrium properties of many body systems [37]. A 
variable range hopping may introduce some degree of disorder thus affecting 
the charge mobility 38, 39J and the thermodynamic functions. 

These issues are hereafter analyzed by means of the path integral formalism 
[40] [4TJ [42] [43] HU [45] which fully accounts for the time retarded e-ph inter- 
actions as a retarded potential naturally emerges in the exact integral action 
|46| . Using the main property of the e-ph coupling model, the Hamiltonian 
linearly dependent on the phonon displacement field, I attack the SSH model 
after introducing a generalized version of the semiclassical model which treats 
only the electrons quantum mechanically. Being valid for any e-ph coupling 
value, the path integral method allows one to derive the partition function and 
the related temperature derivatives without those limitations which affect the 
perturbative studies. The general formalism, both for the ID and 2D system, 
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is described in Section 2 while the path integral approach used to derive the 
full partition function of the interacting system is developed in Section 3. In 
Section 4, 1 outline the main features of the computational method and derive 
the thermodynamical properties for a particle described by a SSH Hamilto- 
nian in a bath of harmonic oscillators. In Section 5, the model is generalized in 
order to include the electron-phonon effects on the oscillator bath: the anhar- 
monic corrections to the total heat capacity are evaluated by a path integral 
cumulant expansion in terms of the source current of the Hamiltonian model. 
Some conclusions are drawn in Section 6. 



2 The Hamiltonian Model 

In a square lattice with isotropic nearest neighbors hopping integral J, the 
SSH Hamiltonian for electrons plus e-ph interactions reads: 

r.s 

+ (Jr,s)y (f r, s+ l f r,s + h.C.) 
(Jr,s)x= --[J + aAu x ] 

(Jr,s)y= --[j + aAUy] 

Au x = u x (r + 1, s) - u x (r, s) 
Au y = u y {r, s + 1) - u y (r, s) , 

(1) 

where a is the electron-phonon coupling, u(r, s) is the dimerization co- 
ordinate indicating the displacement of the monomer group on the (r, s) — 
lattice site, /,t s and / rjS create and destroy electrons (i.e., n band electrons 
in polyacetylene) . The phonon Hamiltonian is given by a set of 2D classical 
harmonic oscillators. The two addenda in Eq. ([T]) deal with one dimensional 
e-ph couplings along the x and y axis respectively, with first neighbors electron 
hopping. Second neighbors hopping processes (with overlap integral J^) may 
be accounted for by adding to the Hamiltonian the term such that 



^ — {J r ,s)x,y{fr+l,s+lfr,s + h.C.') 



(Jr.sky = - - [J (2) + ayJ(Au x ) 2 + (Au y ) 2 



(2) 



The ID SSH Hamiltonian is obtained by Eq. dJ) just dropping the second 
term depending on the y coordinate. The real space Hamiltonian in Eq. ([!]) 
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can be transformed into a time dependent Hamiltonian [47] by introducing the 
electron coordinates: i) (x(r'), y(r')) at the (r,s) lattice site, ii) (x(t), y(r')) 
at the (r + l,s) lattice site and hi) (x(r'), y(r)) at the (r,s + 1) lattice site, 
respectively, r and r' vary on the scale of the inverse temperature /?. The 
spatial e-ph correlations contained in Eq. (fl]) are mapped onto the time axis 
by changing: u x> ( y )(r, s) -> !i Ii(!/ )(t'), u x (r + 1, s) -> ^(t) and u y (r, s + 1) — > 
Uj,(t). Now we set r' = 0, (a;(r'), 2/(r')) = (0,0), {u x (t'),u v (t')) = (0,0). 
Accordingly, Eq. (fTJ) transforms into the time dependent Hamiltonian: 



H(t) = J a (r)(/t( af (r),0)/(0,0) + /».c) 
+ J y (r)(/t(0,y(r))/(0,0) + /i.c. 

^e(t) = -^[^ + a ^( T )] 

J y (r) = --[J + (w s (t)] . 



(3) 



While the ground state of the ID SSH Hamiltonian is twofold degenerate, 
the degree of phase degeneracy is believed to be much higher in 2D [35] as 
many lattice distortion modes contribute to open the gap at the Fermi surface. 
However, as in ID, these phases are connected by localized and nonlinear exci- 
tations, the soliton solutions. Thus, also in 2D both electron hopping between 
solitons [49] and thermal excitation of electrons to band states may take place 
within the model. These features are accounted for by the time dependent 
version of the Hamiltonian. As r varies continuously on the scale and the r- 
dependent displacement fields are continuous variables (whose amplitudes are 
in principle unbound in the path integral), long range hopping processes are 
automatically included in H(t) which is therefore more general than the real 
space SSH Hamiltonian in Eq. |T]) (and Eq. |2])). Thus, by means of the path 
integral formalism we look at the low temperature thermodynamical behavior 
both in ID and 2D searching for those features which may be ascribable to 
some local disorder related to the variable range of the hopping processes. 

The semiclassical nature of the model is evident from Eq. ((3]) in which 
quantum mechanical degrees of freedom interact with the classical variables 
u x( y )( T )- Averaging the electron operators over the ground state we obtain 
the time dependent semiclassical energy per lattice site N which is linear in 
the atomic displacements: 

' J1[T) > = J x (t)P{J,t,x(t)) + J y (j)P(J,T,y(T)) 



N 



P(J,t,v(t)) = i / dkcos[k • v(r)] cosh(ekr)n F (ek) , 



(4) 
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with v(r) = (x(t),0) and v(r) = (0,y(r)) in the first and second term 
respectively, ek = — J J2i= x y cos(ki) is the electron dispersion relation and 
np is the Fermi function. The chemical potential has been pinned to the zero 
energy level. Eq. ^ can be rewritten in a way suitable to the path integral 
approach by defining 

< H{ ^ ] > = V(x(t)) +V(y(r)) +u(r) ■ j(v(r)) 

V(x(r)) = -JP(J,t,x{t)) 
V(y{T)) = -JP(J,r,y(r)) 
j(v(r)) = - a P(j,r,v(r)) 
U ( T ) = (ux{t),u v (t)) . 

(5) 

V(x(t)) and V(y(r)) are the effective terms accounting for the r depen- 
dent electronic hopping while j(v(r)) is interpreted as the external source [3S] 
current for the oscillator field u(r). Averaging the electrons over the ground 
state we neglect the fermion-fermion correlations J50] which lead to effective 
polaron-polaron interactions in non perturbative analysis of the model [51j . 
This approximation however is not expected to affect substantially the follow- 
ing calculations. 



3 The Path Integral Formalism 



Taking a bath of N 2D oscillators, we write the SSH electron path integral, 
C(t) ee (x(r),y(r)), as: 



• exp 



■ exp 



N 



< C(I3) | C(0) >= f[/ Dn t (r) J D((r) 

jyr{ji\T)+V(x{T))+V{y{r))^ 
5>(t).j(v(t)) , 



(6) 



where m is the electron mass, M is the atomic mass and w,; is the oscil- 
lator frequency. As a main feature we notice that the interacting energy is 
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linear in the atomic displacement field. Then, the electronic path integral can 
be derived after integrating out the oscillator degrees of freedom which are 
decoupled along the x and y axis. Accordingly, we get: 



N 



<CC9)IC(o) >= U z 



i=l 



Dx{t) 



exp 
A(x(r)) 



dr{^x 2 (T) + V(x(r))y^A(x(r)) 



.2 N 



E 



AM ^ ftuJi sinh(nwi/3/2) 



drj(x(T)) / dr"co S h(ncj i (\r-r"\-/3/2))j(x(T'% 



2smh(huj i p/2) 



(7) 



Thus, the 2D electron path integral is obtained after squaring the sum 
over one dimensional electron paths. This permits to reduce the computa- 
tional problem which is nonetheless highly time consuming particularly in the 
low temperature limit. Note in fact that the source current j(x(r)) requires 
integration over the 2D Brillouin Zone (BZ) according to Eqs. f4|, ([5|) and 
this occurs for any choice of the electron path coordinates. The quadratic (in 
the coupling a) source action A(x(t)) is time retarded as the particle moving 
through the lattice drags the excitations of the oscillator fields which take a 
time to readjust to the electron motion. When the interaction is sufficiently 
strong the conditions for polaron formation may be fulfilled in the system ac- 
cording to the degree of adiabaticity [S5]. However the present path integral 
description is valid independently of the existence of polarons as it applies also 
to the weak coupling regime. Assuming periodic conditions x(r) = x(t + (3), 
the particle paths can be expanded in Fourier components 



oo 

x(t) = x + 2^R,x n cos(cj„t) — "<sx n sin(u; n T)J 



(8) 



and the open ends integral over the paths J Dx(t) transforms into the 
measure of integration <f> Dx(t). Taking: 



Dx(t) 



dx 



(2Trh 2 /mK B T' 



(1/2) n 



d%tx n d^x n 
(ttH 2 K bT ' i W^) 



(9) 
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we proceed to integrate Eq. {?]) in order to derive the full partition function 
of the system versus temperature. 

Let's point out that, by mapping the electronic hopping motion onto the 
time scale, a continuum version of the interacting Hamiltonian (Eq. ([3])) has 
been de facto introduced. Unlike previous [14] approaches however, our path 
integral method is not constrained to the weak e-ph coupling regime and it 
can be applied to any range of physical parameters. 



4 Computational Method and Thermodynamical Results 

As a preliminar step we determine, for a given path and at a given temper- 
ature: i) the minimum number (iVk) of k— points in the BZ to accurately 
estimate the average interacting energy per lattice site and, ii) the minimum 
number (N T ) of points in the double time integration to get a numerically sta- 
ble source action in Eq. ||7J). The momentum integrations required by Eq. (|4|) 
converge by summing over 1600 and 70 points in the reduced 2D and ID BZ, 
respectively. Moreover, N T = 300 at T = IK. 

Computation of Eqs. jl]) - |7| requires fixing two sets of input parameters. 
The first set contains the physical quantities characterizing the system: the 
bare hopping integral J, the oscillator frequencies u>i and the effective coupling 
X = a 2 h 2 /M (in units meV 3 ). The second set defines the paths for the 
particle motion which mainly contribute to the partition function through: 
the number of pairs (3?x n , 3a;„) in the Fourier expansion of Eq. ([8]), the cutoff 
(A) on the integration range of the expansion coefficients in Eq. © and the 
related number of points (Na) in the measure of integration which ensures 
numerical convergence. 

After introducing a dimensionless path 

N p 

x(r)/a= x + (a n cob^t) + b n sui(o; w t) j , (10) 

with: x~ = x /a, a n = 2^tx n /a and b n = — 2^sx n /a, the functional measure 
of Eq. @ can be rewritten for computational purposes as: 



e2A/a p2A/a p2A/a p2A/a 

dai / dbi ■ ■ ■ da,N p / db~N p 

-2A/a J-2A/a J -2A/a J -2A/a 



(11) 



In the following we take the lattice constant a = lA. As a criterion to set 
the cutoff A on the integration range, we notice that the functional measure 
normalizes the kinetic term in Eq. ([7]): 
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Dx(r)exp —— / cItx 2 (t) 



1 



(12) 



and this condition holds for any number of pairs N p truncating the Fourier 
expansion in Eq. ©. Then, taking N p = 1, the left hand side of Eq. ([TT 
transforms as: 




T Jo dT±2[T) . 



U I 2 
dyexp{-y 2 ) 



ijiKbT 



Using the series representation 



/ dyexp(-y 2 ) = 
Jo fc=0 



(-l) fc C/ 



2fe+l 



fc!(2fc+ 1) 



(13) 



(14) 



one determines U (after setting the series cutoff k max which ensures con- 
vergence) by fitting the Poisson integral value s/tt/2. 

Thus, we find that the cutoff A can be expressed in terms of the thermal 
wavelength A as A ~ A/V27T 3 hence, it scales versus temperature as A oc 
1/y/T. This means physically that, at low temperatures, A is large since many 
paths are required to yield the correct normalization. For example, at T = IK, 
we get A ~ 284A Numerical investigation of Eq. (0 shows however that a 
much shorter cutoff suffices to guarantee convergence in the path integral, 
while the cutoff temperature dependence implied by Eq. ([12"]) holds also in 
the computation of the interacting partition function. The thermodynamical 
results hereafter presented have been obtained by taking A ~ A/(10\/27r 3 ). 
Summing in the ID system over Na ~ 20/VT points for each integration range 
and taking N p = 2, we are then evaluating the contribution of (Na + \^ 2N p+ 1 
paths (the integer part of JV^ is obviously selected at any temperature). Thus, 
at T = IK and in ID, we are considering ~ 4 • 10 6 different paths for the 
particle motion while, at T = 100A" the number of paths drops to 243 [54]. In 
the 2D problem, Na ~ 35/VT points for each Fourier coefficient are required 
[55] , Computation of the second order derivatives of the free energy in the 
range T S [1, 301]^, with a spacing of 3K, takes 55hours and 15 minutes on a 
Pentium 4. Note that larger N p in the path Fourier expansion would further 
increase the computing time without introducing any substantial improvement 
in the thermodynamical output of our calculation. 

Although the history of the SSH model is mainly related to wide band poly- 
mers, we take here a narrow band system (J = lQOmeV) [28] with the caveat 
that electron-electron correlations may become relevant in narrow bands. 
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Free energy and heat capacity have been first computed up to room tem- 
perature, both in ID and 2D, assuming a bath of N — 10 low energy os- 
cillators separated by 2meV: %u>\ = 2meV, Tk^io = 20meV. The lowest 
energy oscillator yields the largest contribution to the phonon partition func- 
tion mainly in the low temperature regime while the ui±o oscillator essentially 
sets the phonon energy scale which determines the size of the e-ph coupling. 
A larger number N of oscillators in the aforegiven range would not signifi- 
cantly modify the calculation whereas lower Wi values would yield a larger 
contribution to the phonon partition function mainly at low T. 

In the discrete SSH model, the value a = Aa 2 / (ttkJ) ~ 1, marks the 
crossover between weak and strong e-ph coupling, with k being the effective 
spring constant. In our continuum and semiclassical model the effective cou- 
pling is X- Although in principle, discrete and continuum models may feature 
non coincident crossover parameters, we assume that the relation between a 
and J obtained by the discrete model crossover condition still holds in our 
model. Hence, at the crossover we get: \ c ~ 7rJ7i wf /64. This means that, 
in Figures 1-3, the crossover value is set at \c ~ 2000meV 3 . 

In Fig.l, a comparison between the ID and the 2D free energies is presented 
for two values of %, one lying in the weak and one in the strong e-ph coupling 
regime. The oscillator free energies F p h are plotted separately while the free 
energies arising from the total action in Eq. (|7jl. shortly termed F source , results 
from the competition between the free path action (kinetic term plus hopping 
potential in the exponential integrand) and the source action depending on the 
e-ph coupling. While the former enhances the free energy, the latter becomes 
dominant at increasing temperatures thus reducing the total free energy. In 
general, the 2D free energies have a larger gradient (versus temperature) than 
the corresponding ID terms. The F p h lie above the F source both in ID and 2D 
because of the choice of the Juvf. in general one may expect a crossing point 
between F p h and F source with temperature location depending on the value 
of X- Note that the F source plots have a positive temperature derivative in the 
low temperature regime and this feature is more pronounced in 2D. In fact, 
at low T, the source action is dominated by the hopping potential V(x(r)) 
while, at increasing T, the e-ph effects become progressively more important 
(as the bifurcation between the xi and X2 curves shows) and the F source get 
a negative derivative. In 2D, the weight of the VCx(t)) term is larger be- 
cause there is a higher hopping probability. This physical expectation is taken 
into account by the path integral method. At any temperature, we monitor 
the ensamble of relevant particle paths over which the hopping potential is 
evaluated. For a selected set of Fourier components in Eq. (fl"U|) the hopping 
decreases by lowering T but its value is still significant at low T. Considered 
that: i) the total action is obtained after a dr integration of V[x(t)j and ii) 
the dr integration range is larger at lower temperatures, we explain why the 
overall hopping potential contribution to the total action is responsible for 
the anomalous free energy behavior at low T. 
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In Fig. 2, the heat capacity contributions due to the oscillators (C p h) and 
electrons plus e-ph coupling (C sou ) are reported on. The values are normalized 
over N. the previously described summation over a large number of paths 
turns out to be essential to recover the correct thermodynamical behavior 
in the zero temperature limit. The dimensionality effects are seen to be large 
and, for a given dimensionality, the role of the e-ph interactions is magnified at 
increasing T. The total heat capacity (C p h + C SO u) overT ratios are plotted in 
Fig. 3 in the low T regime to emphasize the presence of an anomalous upturn 
which appears at T < 10K in ID and T < 20K in 2D. This feature in the heat 
capacity linear coefficient is ultimately related to the sizable effective hopping 
integral term V(a;(r)). The strength of the e-ph coupling has a minor role in 
the low T limit although it determines the shape of the anomaly versus T. 

Let's focus now on the ID system and consider the effect of the oscillators 
bath on the thermodynamical properties: in Figures 4-6 the ten phonon en- 
ergies are: Uui\ — 22meV 7 huJio — AOmeV. Accordingly the crossover is set 
at Xc ~ SOOOmeV 3 and three plots out of five lie in the strong e-ph coupling 
regime. As shown in Fig. 4, large x values are required to get strongly de- 
creasing free energies versus temperature while the x = SOOOmeV 3, curve now 
hardly intersects the phonon free energy at room temperature. Fig. 5 shows 
the rapid growth of the source heat capacity versus temperature at strong 
couplings whereas the presence of the low T upturn in the total heat capacity 
over T ratio is confirmed in Fig. 6. Note that, due to the enhanced oscillators 
energies, the phonon heat capacity saturates here at T ~ 400K (in Fig. 2, for 
the ID case, T ~ 200K). 

5 Electron-Phonon Anharmonicity 

So far we have considered a bath of harmonic phonons. Now we face the 
following question: what is the effect of the particle-phonon interaction on the 
phonon subsystem? 

In general, the phonon partition function perturbed by a source current 
j(r) can be expanded in anharmonic series as: 

Z ph \j(r)]^ Z ph h +£(_!)' <C >. (T) 
^ i=i 

(15) 

where the cumulant terms < C l >j( r ) are expectation values of powers of 
correlation functions of the perturbing potential. The averages are meant over 
the ensemble of the harmonic oscillators whose partition function is Z p h- 
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Fig. 1. Phonon (F p h) and Source Term (F aou ) contributions to the ID and 2D 
free energies for two values of the effective coupling x- Xi ~ 1440meV 3 (weak e-ph 
coupling) \2 = 2560rneV 3 (strong e-ph coupling). A bath of ten phonon oscillators 
is considered, the largest phonon energy being huio = 20meV. 
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Fig. 2. Phonon and Source Term contributions (normalized over the number of 
oscillators) to the ID and 2D heat capacities for the same parameters as in Fig.l. 
The oscillator heat capacities are also plotted. 
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Fig. 6. Total heat capacity over temperature for the same parameters as in Fig. 4. 
The phonon heat capacity is also plotted. 
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5.1 The Holstein-type Current 

First we consider the general problem of an electron path linearly coupled to 
a single oscillator with energy u and displacement u(r), through the current 
jx(j) = —ax(r). This type of current models the Holstein interaction [56] . In 
this case odd k cumulant terms vanish and the lowest order even k cumulants 
can be straightforwardly derived [S7] . To obtain a closed analytical expression 
for the cumulants to any order we approximate the electron path by its r 
averaged value: < x(t) >= drx(r) = xo/a and expand the oscillator 

path in Np Fourier components: 



u(t) = u + ^2 ^(jfcun cos(cj„t) — $sx n sin(o; n T)^ 

n=l 

uj n = 2irn/f3 . (16) 
Next we choose the measure of integration 



2N F 



(2tt ■ -2N f tt 



DU[T) = 1 2 ) ^Aif- 1 ) ' ^ 



X 



Yl / Mu n / d%U n , (17) 



being Am = \Jirti 2 [3/M. Such a measure normalizes the kinetic term in 
the oscillator field action 



M 



Du(r)exp - — / <1tu 2 (t) = 1 . (18) 



o 



Then, using Eqs. (p~5]) - (fTTj) . we obtain for the k — th cumulant 



_ v -i (<* fl /?A M ) fc (fc-l)!! fff (2mr) 2 



ph k\ir k / 2 (u>l3) k+1 J-JL (2 n7 r) 2 + (w/3) 2 
a_R = ax /a . (19) 

Let's set Xo/a = 0.1 in the following calculations thus reducing the effec- 
tive coupling an by one order of magnitude with respect to the bare value. 
However, the trend shown by the results hereafter presented does not depend 
on this choice since xq /a and a can be varied independently. As the cumulants 
should be stable against the number of Fourier components in the oscillator 
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path expansion, using Eq. (|19p we set the minimum Np through the condi- 
tion 2NpTT ^> uf3. The thermodynamics of the anharmonic oscillator can be 
computed by the cumulant corrections to the harmonic phonon free energy: 



F(*'(T) 



1 

1 



In 



k 

£ 

i=i 



< C 21 >N F 



(20) 



To proceed one needs a criterion to find the temperature dependent cutoff 
k* in the cumulant series. We feel that, in the low T limit, the third law of ther- 
modynamics may offer the suitable constraint to determine k* . Then, given 
a and uj, the program searches for the cumulant order such that the phonon 
heat capacity and the entropy tend to zero in the zero temperature limit. At 
any finite temperature T, the constant volume heat capacity is computed as 



C« ( T ) = _ F W (rp + 2A) _ 2F (k) ( T + A)+ p(k) ( T ) 

A being the incremental step and k* is determined as the minimum value 
for which the heat capacity converges with an accuracy of 10~ 4 . Figs. 7(a) and 
7(b) show phonon heat capacity and free energy respectively in the case of 
a low energy oscillator for an intermediate value of e-ph coupling. Harmonic 
functions, anharmonic functions with second order cumulant and anharmonic 
functions with k* corrections are reported on in each figure. The second or- 
der cumulant is clearly inadequate to account for the low temperature trend 
yielding a negative phonon heat capacity below ~ AOK while at high T the 
second order cumulant contribution tends to vanish. Instead, the inclusion of 
k* terms in Eqs. ([20)) . (|2Tj) leads to the correct zero temperature limit although 
there is no visible anharmonic effect on the phonon heat capacity throughout 

(k*) 

the whole temperature range being C v perfectly superimposed on the har- 
monic Gy. Note in Fig. 7(b) that the k* corrections simply shift downwards 
the free energy without changing its slope versus temperature. By increasing 
an, the low T range with wrong (negative) C v broadens whereas the k* con- 
tributions permit to fulfill the zero temperature constraint and substantially 
lower the phonon free energy. Thus, for the particular choice of constant (in r) 
source current we find that the e-ph anharmonicity renormalizes the phonon 
partition function although no change occurs in the thermodynamical behav- 
ior of the free energy derivatives. Anharmonicity is essential to stabilize the 
system but it leaves no trace in the heat capacity [58]. Figure 8(a) displays 
the k* temperature dependence for three choices of e-ph coupling in the case 
of a low energy oscillator: while, at high T, the number of required cumulants 
ranges between six and ten according to the coupling, k* strongly grows at low 
temperatures reaching the value 100 at T — IK for ccr = GOmeVA -1 . The 
k* versus olr behavior is depicted in Fig. 8(b) for three selected temperatures: 
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at low T the cutoff strongly varies with the strength of the coupling while, 
by enhancing T, the number of cumulant terms in the series is smaller and 
becomes much less dependent on a^. 



5.2 The SSH-type Current 

Next we turn to the computation of the equilibrium thermodynamics of the 
phonon subsystem perturbed by the source current of the semiclassical SSH 
model described in Sections 2 and 3. Assuming that the electron particle 
path interacts with each of the N oscillators through the coupling a (taken 
independent of i), we write the k th cumulant term as 



(22) 

where j(r) is given in Eq. (0. We take here the ID system. Since the oscil- 
lators are in fact decoupled in our model (and anharmonic effects mediated by 
the electron particle path are neglected) the behavior of the cumulant terms 
< C k >j( T ) can be studied by selecting a single oscillator having energy ui and 
displacement u(t). 

As the electron propagator depends on the bare hopping integral we set as 
before J = lOOmeV. Any electron path yields in principle a different cumulant 
contribution. Numerical investigation shows however that convergent k— order 
cumulants are achieved by taking N p = 2 Fourier components in the electron 
path expansion and summing over ~ 5 2Np+1 electron paths. 

As in the case of Sect. 5.1, we truncate the cumulant series by invok- 
ing the third law of thermodynamics to determine the cutoff k* in the low 
temperature limit and by searching for numerical convergence on the first 
and second free energy derivatives at any finite temperature. Again, we can 
start our analysis from Eq. (|2"0|) after checking that odd k cumulants yield 
vanishing contributions. Now however the picture of the anharmonic effects 
changes drastically. The e-ph coupling strongly modifies the shape of the heat 
capacity and free energy plots with respect to the harmonic result as it is 
seen in Figs. 9(a) and 9(b) respectively. The heat capacity versus temperature 
curves show a peculiar peak above a threshold value a ~ lOmeV A -1 which 
clearly varies according to the energy of the harmonic oscillator. Here we set 
uj = 20meV to emphasize the size of the anharmonic effects on a low en- 
ergy oscillator. By enhancing a the height of the peak grows and the bulk of 
the anharmonic effects on the heat capacity is shifted towards lower T. At 



n k r p 

<C*>i(r)= z pl !Hf Du t (r)-J[ / dnMn)j{n) 



x exp 



ff 3 N M- 
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a ~ 60meV 'A -1 the crossover temperature is around 100K. Note that the size 
of the anharmonic enhancement is ~ 10 times the value of the harmonic oscil- 
lator heat capacity at T = 100A'. It is worth noting that previous numerical 
studies of a classical one dimensional anharmonic model undergoing a Peierls 
instability [59] also found a specific heat peak as a signature of anharmonicity. 
However such a large anharmonic effect on the phonon subsystem is partly 
covered in the total heat capacity by the source action A(j(r)) and mainly by 
the hopping potential V(x(t)) contributions analysed in the previous Section. 

Taken for instance a bath of ten low energy oscillators with lu = 20meV, 
setting a ~ QOmeVA -1 which implies an effective coupling \ ~ 700meV 3 (last 
of Eqs. ([5])) we get a source heat capacity a factor two larger than the harmonic 
phonon heat capacity at temperatures of order 100K . Thus the anharmonic 
peak, although substantially smeared by the electronic contributions to the 
total heat capacity, should still appear in systems with low energy phonons 
and sizeable e-ph coupling to which the SSH Hamiltonian applies. Let's focus 
on this point. 

The total heat capacity is given by the phonon contribution plus a source 
heat capacity which includes both the electronic contribution (related to the 
electron hopping integral) and the contribution due to the source action (the 
latter being oc a 2 ). Fig. 10(a) shows the comparison between the anharmonic 
phonon heat capacity (C v nh ) and the source heat capacity, here termed Cy ~ p 
to emphasize the dependence both on the electronic and on the e-ph coupling 
terms. C v ~ p is computed as described above setting a — 21.74meVA~ 1 [55] . 
Also the total heat capacity (C£?* = C v nh + Cy~ p ) is shown in Fig.lO(a). 
At low temperatures, CfT p yields the largest effect mainly due to the elec- 
tronic hopping while at high T, Cy~ p prevails as the source action becomes 
dominant. In the intermediate range (T G [90,210]) the anharmonic phonons 
provide the highest contribution although their characteristic peak is sub- 
stantially smeared in the total heat capacity by the source term background. 
Fig. 10(b) compares Cy 1 and C v ~ p for two increasing values of a: while the 
anharmonic peak shifts downwards (along the T axis) by enhancing a, C v nh 
remains larger than C v ~ p in a temperature range which progressively shrinks 
due to the strong dependence of the source action on the strength of the e-ph 
coupling. Finally, we observe that the low temperature upturn displayed in 
the total heat capacity over T ratio discussed above is not affected by the 
inclusion of phonon anharmonic effects which tend to become negligible at 
low temperatures. 

6 Conclusions 

Mapping the real space Su-Schrieffer-Heeger model onto the time scale I have 
developed a semiclassical version of the interacting model Hamiltonian in 
one and two dimensions suitable to be attacked by path integrals methods. 
The acoustical phonons of the standard SSH model have been replaced by a 
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set of oscillators providing a bath for the electron interacting with the dis- 
placements field. Time retarded interactions are naturally introduced in the 
formalism through the source action A(x(t)) which depends quadratically on 
the bare e-ph coupling strength a. Via calculation of the electronic motion 
path integral, the partition function can be derived in principle for any value 
of a thus avoiding those limitations on the e-ph coupling range which bur- 
den any perturbative method. Particular attention has been paid to establish 
a reliable and general procedure which allows one to determine those input 
parameters intrinsic to the path integral formalism. It turns out that a large 
number of paths is required to carry out low temperature calculations which 
therefore become highly time consuming. The physical parameters have been 
specified to a narrow band system and the behavior of some thermodynamical 
properties, free energy and heat capacity, has been analysed for some values 
of the effective coupling strength lying both in the weak and in the strong 
coupling regime. We find, both in ID and 2D, a peculiar upturn in the low 
temperature plots of the heat capacity over temperature ratio indicating that 
a glass-like behavior can arise in the linear chain as a consequence of a time 
dependent electronic hopping with variable range. 

According to our integration method (Eq. (|T|)). at any temperature, a 
specific set of Fourier coefficients defines the ensamble of relevant particle 
paths over which the hopping potential y(x(r)) is evaluated. This ensamble 
is therefore T dependent. However, given a single set of path parameters one 
can monitor the V(x(t)) behavior versus T. I find that the hopping decreases 
(as expected) by lowering T but its value remains appreciable also at low 
temperatures 20K in 2D and ■< 10K in ID). Since the dr integration range 
is larger at lower temperatures, the overall hopping potential contribution 
to the total action is relevant also at low T. It is precisely this property 
which is responsible for the anomalous upturn in the heat capacity linear 
coefficient. Further investigation also reveals that the upturn persists both in 
the extremely narrow (J ~ lOmeV) and in the wide band (J ~ leV) regimes. 
Moreover, the upturn is not modified by the inclusion of electron-phonon 
anharmonicity in the phonon subsystem. 

The presented computational method accounts for a variable range hop- 
ping on the r scale which corresponds physically to introduce some degree of 
disorder along the linear chain. This feature makes my model more general 
than the standard SSH Hamiltonian (Eq. fl}) with only real space nearest 
neighbor hops. While hopping type mechanisms have been suggested [49] to 
explain the striking conducting properties of doped polyacetylene at low tem- 
peratures I am not aware of any other direct calculation of the specific heat in 
the SSH model. Since the latter quantity directly probes the density of states 
and integrating over T the specific heat over T ratio one can have access to 
the experimental entropy, this method may provide a new approach to anal- 
yse the transition to a disordered state which indeed exists in polymers. In 
this connection it is also worth noting that the low T upturn in the specific 
heat over T ratio is a peculiar property of glasses |60|, 16 lj in which tunneling 
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states for atoms (or group of atoms) provide a non magnetic internal degree 
of freedom in the potential structure [551 [53] . 
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Fig. 7. (a) ID Phonon heat capacity and (b) ID Phonon free energy calculated 
in i) the harmonic model, ii) anharmonic model with second order cumulant, iii) 
anharmonic model with k* cumulants (see text), an is the effective e-ph coupling in 
units meV A 1 and w is the phonon energy. 
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Fig. 8. (a) Number of cumulants required to obtain a convergent phonon heat 
capacity at any temperature for different choices of e-ph couplings, (b) Number of 
cumulants yielding a convergent phonon heat capacity at any e-ph coupling for three 
selected temperatures. 
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Fig. 9. ID Anharmonic (a) phonon heat capacity and (b) free energy versus tem- 
perature for eight values of e-ph coupling. The harmonic plots are also reported on. 
A low energy oscillator is assumed. 
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Fig. 10. (a) Total Heat Capacity versus temperature in the ID Su-Schrieffer-Heeger 
model. The contributions due to anharmonic phonons (C v nh ) and electrons plus 
electron-phonon interactions (Cy~ p ) are plotted separately. The largest a of Fig.l 
is assumed, (b) C v nh and Cy V for two values of a in units meVA^ 1 . ui = 20meV. 
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